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System analysis and design of any entry system must balance the level fidelity for each 
discipline against the project timeline. One way to inject high fidelity analysis earlier in the 
design effort is to develop surrogate models for the high-fidelity disciplines. Surrogate 
models for the Thermal Protection System (TPS) are formulated as Mass Estimating 
Relationships (MERs). The TPS MERs are presented that predict the amount of TPS 
necessary for safe Earth entry for blunt-body spacecraft using simple correlations that 
closely match estimates from NASA’s high-fidelity ablation modeling tool, the Fully Implicit 
Ablation and Thermal Analysis Program (FIAT). These MERs provide a first order estimate 
for rapid feasibility studies. There are 840 different trajectories considered in this study, and 
each TPS MER has a peak heating limit. MERs for the vehicle forebody include the ablators 
Phenolic Impregnated Carbon Ablator (PICA) and Carbon Phenolic atop Advanced 
Carbon-Carbon. For the aftbody, the materials are Silicone Impregnated Reusable Ceramic 
Ablator (SIRCA), Acusil II, SLA-561V, and LI-900. The MERs are accurate to within 14% 
(at one standard deviation) of FIAT prediction, and the most any MER under predicts FIAT 
TPS thickness is 18.7%. This work focuses on the development of these MERs, the resulting 
equations, model limitations, and model accuracy. 


cw 

= cold wall 

HL 

= total heat load, J/cm 2 

MER 

= mass estimating relationship 

P 

= surface pressure, Pa 

SD 

= standard deviation 

TH 

= TPS thickness, cm 

TPS 

= thermal protection system 

V 

= entry velocity, km/s 

Voo 

= free stream velocity, km/s 

Poo 

= free stream density, kg/m 3 


Nomenclature 


I. Introduction 

A multidisciplinary, integrated tool called the “Multi Mission System Analysis for Planetary Entry Descent and 
Landing,” also known as M-SAPE 1 , is being developed as part of the Entry Vehicle Technology project under 
NASA’s In-Space Propulsion Technology program. The primary purpose of M-SAPE is to perform system analysis 
and design for an Earth entry vehicle suitable for sample return missions. It includes geometry, mass sizing, impact 
analysis, structural analysis, flight mechanics, TPS, and a web portal for user access. It also provides a platform by 
which technologies and design elements can be evaluated rapidly prior to any costly investment commitment. An 
important feature of for M-SAPE is the ability to perform system analysis and trade studies in minutes and hours 
(not weeks or months). Part of M-SAPE’s application requires the development of parametric mass estimating 
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relationships (MERs) to determine a vehicle’s required Thermal Protection System (TPS). For this analysis, the heat 
shield is made of a uniform thickness TPS, and the resulting MERs determine the pre-flight mass of the TPS. 

The analysis and design of a vehicle’s safe Entry, Descent, and Landing (EDL) is very multidisciplinary in 
nature, requiring the application of Guidance, Navigation, and Control (GNC), aerodynamics, aerothermo dynamics, 
and material response. A schematic diagram of the process is given in Fig. 1. For a typical re-entry problem, 
designers would first choose a payload, vehicle geometry, and subsequent vehicle mass. A trajectory is then 
determined by the use of GNC, with the trajectory uniquely identified by its entry velocity, entry flight path angle, 
and ballistic coefficient. The vehicle’s aero thermal environment is then found by the use of computational 
aero thermodynamics. An aero thermal analysis will include convective and radiative heating, surface pressure, and 
recovery enthalpy. Once the entry environment is known, the TPS material response is modeled to determine the 
material thickness necessary to keep the bond line temperature below a specified value. It is also important to know 
the amount of surface recession, if any. 

The traditional approach for this coupled problem would be to use a high fidelity computational fluid dynamics 
code such as the Data Parallel Line Relaxation 2 (DPLR) code or the Langley Aeroheating Upwind Relaxation 
Algorithm 3 (LAURA) for the aero thermal component. For TPS response one could use the Fully Implicit Ablation 
and Thermal Response 4 (FIAT) code, the Charring Material Thermal Response and Ablation Program 5 (CMA), the 
Charring Ablating Thermal Protection Implicit System Solver 6 (CHAR), or the Standard Ablation Program 7 
(STAB). However, this coupled approach usually has a very slow turnaround time and is highly dependent upon 
analyst availability. 

The MERs developed in this work are to be used as alternates to high fidelity TPS material response codes such 
as FIAT or CHAR. The goal of the current work is to use these MERs in M-SAPE to screen a proposed trajectory to 
see if the vehicle’s resulting TPS mass is too much or if a trajectories’ heat rate exceeds a TPS limit. 



Figure 1. Flow chart of the EDL design process 

For the vehicle forebody stagnation point, M-SAPE uses the Sutton - Graves 8 correlation for convective heating 
and the Tauber - Sutton 9 correlation for radiative heating. To date, however, no correlations based on high-fidelity 
FIAT modeling are known. As a result, the current work is to develop FIAT-based MERs that match FIAT 
prediction as closely as possible. Six MERs are developed. For the vehicle forebody the ablative materials are 
Phenolic Impregnated Carbon Ablator 10 (PICA) and Carbon Phenolic 11 atop Advanced Carbon-Carbon. For the 
vehicle aftbody the materials are Silicone Impregnated Reusable Ceramic Ablator 12 (SIRCA), Acusil 13 II, SLA- 
561V, and LI-900. As will be shown, the MERs are accurate to FIAT prediction within 14% at one standard 
deviation. For these MERs, no margins have been added to the TPS thickness. 

Using these MERs, M-SAPE can now perform rapid trade studies involving entry velocity, ballistic coefficient, 
vehicle geometry, entry flight path angle, etc. for required TPS thickness and mass. Using these MERs have reduced 
design turnaround times for a possible Earth entry configuration from weeks to hours. 

II. MER Development 
A. Concern Regarding the Application of MERs 

The MERs are statistical correlations developed to predict FIAT output. Each MER has with it a listed accuracy 
to the FIAT code’s prediction, which is the standard deviation of the MER/FIAT thickness ratio over all sample 
trajectories. 

The MER for each material has a maximum allowable instantaneous heat flux (combined convective and 
radiative). M-SAPE has correctly implemented this restriction into its code. If a trajectory’s peak heating (found 
from the Sutton-Graves and Tauber-Sutton correlations) exceeded a TPS material’s allowed value (see Table 4), 
then the trajectory was discarded and not included in the analysis. 

It is emphasized that the MER thickness is not the manufacturing limit, and that substantial margin may be 
added. For example, Stardust flew with a PICA thickness 14 of 5.816 cm (2.29 inch), which is much greater than its 
unmargined thickness. 
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Finally, as for any statistical analysis there exist some extreme trajectories for which the ratio of MER/FIAT TPS 
thickness prediction can well exceed the listed MER accuracy at one SD. For the present work, the largest MER 
under prediction of FIAT is 18.7%. It is for this reason that full datasets are presented showing the curve fit and 
MER/FIAT “goodness-of-fif ’ data. M-SAPE uses these MERs as a rough approximation to determine TPS thickness 
for flight trajectories of interest, but a true high-fidelity analysis is always a requirement for a proposed mission. 

B. Flight Trajectory Parameters 

Found in Table 1 is information on the flight trajectory parameters used for this study. The range covers most of 
sample return applications from the Moon, Mars, and comets that are close to Earth orbit. The level accuracy reflects 
the needs for a conceptual level design. 


Table 1. Flight trajectories considered for the MERs 


Flight Trajectory Parameter 

Range of Values 

Resolution 

Entry Velocity [km/s] 

10-16 

1 

Entry Flight Path Angle [abs. deg.] 

5-25 

5 

Ballistic Coefficient [kg/m 2 ] 

41.95- 128.74 

15.5 (max) 

Total number of trajectories 

840 

- 


Here, resolution is defined as a parameter’s smallest step-size 

C. Vehicle Geometry Parameters 

The vehicle is a 60° sphere-cone shape. Details of the vehicle geometry are given in Table 2. 

Table 2. Vehicle Geometry 



min 

max 

Nose Radius [m] 

0.352 

0.352 

Vehicle Diameter [m] 

0.75 

2.25 

Vehicle Mass [kg] 

16 

159 

Payload Mass [kg] 

5 

30 


D. FIAT modeling constraints 

FIAT analysis of each trajectory uses the following constraints: 

• The maximum temperature at the bottom face of the top material (bond line) is 250°C. 

• The back face of the material stackup is adiabatic. 

• The surrounding environment is at 21.3°C for radiation calculations from the spacecraft surface. 

• ID planar geometry 

• No margins are added to the thickness 

• FIAT v3.0 

It should be noted that FIAT is a ID code and is most applicable for heat shield regions that do not undergo rapid 
shape changes leading to large variations in aeroheating with spatial location. Examples of these regions include the 
stagnation point, along the flank, or any other acreage location. For regions that do change shape quickly, such as at 
the shoulder, a material response code like TITAN, 15 3dFIAT, or CHAR is more appropriate because it includes 2D 
or 3D effects. In addition, PICA’s heat conduction is orthotropic, which also necessitates the use of multi- 
dimensional codes along regions of the heat shield that change shape quickly. Shown in Table 3 are the ranges of 
sample trajectory heating rates, heat loads, and surface pressures used in the FIAT analysis of the 840 sample flight 
trajectories. The data shown in this table are prior to MER development. 

Table 3. Range of surf ace heating, heat load, and surface pressure. 



Forebody 

Aftbody 

Maximum heat flux [W/cm 2 ] 

151 -3767 

2-58 

Heat Load [J/cm 2 ] 

3855-34453 

59-531 

Maximum pressure, atm 

0.03-3.18 

0.003 - 1.64 


E. Methodology 
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1. MER Formulation 

The variables considered to create the MERs included peak heat flux, peak surface pressure, heat load, ballistic 
coefficient, entry velocity, and entry flight path angle. To determine which of these variables to use, sensitivity 
studies were conducted by plotting the value of each variable against required TPS thickness for all trajectories. 
Viewing these scatter plots and calculating correlation coefficients determined if any correlation or sensitivity 
existed. For this work, heat load and entry velocity were chosen for the MER correlations. We have selected a ratio 
of heat load over square of velocity as the correlation parameter that represents the ratio of entry thermal energy 
over entry kinetic energy. The resulting MER accuracy is sufficient to meet M-S APE’s needs. 

2. Forebody Calculations 

The form of the Sutton-Graves relation used by M-SAPE is a cold-wall (CW) convective heat flux. 
Consequently, the convective heat load used in the MERs is CW. For the radiative heat load, the Tauber- Sutton 
relation was used to estimate the radiative heating. Stagnation pressure was found using the momentum equation: 

P = PooV£ (1) 

When FIAT ran through all 840 trajectories, 123 of them were so mild as to produce no recession or ablation for 
PICA and Carbon Phenolic. For such mild environments, the proper heat shield material is an insulator rather than 
an ablator. As a result, these mild trajectories were not used in the forebody MERs, and the MER minimum 
thickness was chosen from the trajectories that cause any ablation. This minimum thickness is enough material to 
keep the bond line temperature below 250°C for all trajectories. 

3. Aftbody Calculations 

To estimate the aero thermal environment on the aftbody, convective heating was found by taking 5% of the 
forebody stagnation point heating. This value is consistent with a study showing Viking 16 aftbody (base) heating to 
be about 5% of Sutton-Graves forebody heating. On the aftbody, radiative heating was ignored. Aftbody surface 
pressures are suggested to be 30% of forebody based on Earth re-entry flight test data 17 using slender cones with 
blunted nose tips and a nearly flat base. For this work, a surface pressure of 50% of the forebody stagnation point 
value was used as a conservative estimate. 

F. Summary of Results 

The results are summarized in Table 4. A discussion of these results, how they were obtained, and their accuracy 
is given in Section III. 

T able 4. Summary of results. 



Forebody 

Aftbody | 

PICA 

CP/ACC 

SIRCA 

Acusil II 

LI-900 

SLA-561V 

Max allowable 
heat flux, CW 
[W/cm 2 ] 

1200 

30000 

100 

100 

75 

100 

Recession [cm] 

0.60-1.15 

0.0002-0.141 

none 

none 

none 

none 

Accuracy to 
FIAT [one SD] 

6.3% 

7.3% 

8.5% 

7.6% 

14.0% 

8.5% 

Largest under 
prediction [% of 
FIAT] 

11.7 

16.6 

15.6 

15.1 

18.7 

15.7 

Minimum 
thickness [cm] 

3.27 

2.27 

0.518 

0.614 

0.686 

0.454 

Thickness [cm] 


/ m \ 0 ' 2102 

U959(f) 

o, 281 (tr 

( rrr \ 0-5697 

q cw <= 10 W/cm 2 

q cw > 10 W/cm 2 
-0.0306 +0.5896^j + 0.6739 

0.0064^j +0.0961^1 + 0.3322 


III. Forebody MER Development: PICA and Carbon Phenolic 

7. PICA 

Phenolic Impregnated Carbon Ablator (PICA) is a low-density, high performance ablator that uses a porous rigid 
carbon-fiber substrate impregnated with a dispersed phenolic polymer. It has a nominal density of 0.274 g/cm 3 (17.1 
lb/ft 3 ). 18 PICA has successfully flown as the main forebody heat shield material on the Stardust Sample Return 
Capsule (SRC) 19 and on the Mars Science Laboratory (MSL) capsule. 20 
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The material stackup for this MER consists of only PICA, and the FIAT PICA material response model was 
version 3.3. The MER is given in Eq. 2, and it has a standard deviation of 6.3% to FIAT prediction. Specifics of the 
correlation model are given in Table 5. Figure 2 shows the MER fit of the data, and a scatter plot of Goodness of Fit 
(GoF) with FIAT-predicted TPS thickness is given in Fig. 3. GoF is defined as the ratio of MER to FIAT prediction. 

TH = 1.8686 (Jf J (2) 


Table 5. PICA MER Details 


Variable 

Values 

Maximum allowable heat flux, W/cm 2 , CW 

1200 

Recession, cm 

0.134-1.153 

Accuracy to FIAT at one SD 

6.3% 

Fargest under prediction of FIAT (% of FIAT) 

11.7 

Number of trajectories with no recession 

123 

Trajectories used for correlation 

419 

Minimum thickness, cm 

3.27 



Figure 2. PICA MER banded by 1 standard deviation 



45 th AIAA Thermophysics Conference, Dallas Texas 

June 22-26, 2015 



Figure 3. PICA MER Goodness of Fit 


2. Carbon Phenolic Atop Advanced Carbon-Carbon 6 

Carbon Phenolic is a high-density ablator capable of withstanding extremely high heat flux environments. The 
maximum allowable heat flux for this material is an amazing 30 kW/cm 2 , which it approximately experienced on the 
Galileo Jupiter Probe planetary entry. 21 Carbon phenolic has also been used on the four Pioneer Venus probes. 22 
There are different methods to fabricate this material, and the one used for this work was specified as FM-5055. 23 
The virgin material density 24 is 1.43 g/cm 3 (89.6 lb/ft 3 ). 

For this MER, the material stackup is given in Table 6, with ACC6 as the carrier structure. 

Table 6. Material sta ckup 


Material 

Thickness, cm 

Carbon Phenolic 

variable 

HT-424 (adhesive) 

0.0381 

Advanced Carbon-Carbon (ACC) version 6 

0.250 


This MER is given in Eq. 3, and it has a standard deviation of 7.3% to FIAT prediction. Specifics of the MER 
are given in Table 7. A scatter plot showing the fit to the data is given in Fig. 4, and Fig. 5 shows a scatter plot of 
GoF with FIAT-predicted TPS thickness. 


TH = 1.1959 



0.2102 


( 3 ) 


Table 7. Carbon P henolic atop ACC MER Details 


Variable 

Values 

Maximum allowable CW heat flux, W/cm 2 

30000 

Recession [cm] 

0.0002-0.1416 

Accuracy to FIAT at one SD 

7.3% 

Largest under prediction of FIAT (% of FIAT) 

16.6 

Number of FIAT non-convergent trajectories 

16 

Trajectories with no recession 

123 

Trajectories used for correlation 

701 

Minimum thickness, cm 

2.266 



Figure 4. Carbon Phenolic over ACC6 MER banded by 1 standard deviation. 
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Figure 5. Carbon Phenolic over ACC6 MER GoF. 


IV. Vehicle Backshell MERs 

1 . SIRCA 

Silicone Impregnated Reusable Ceramic Ablator (SIRCA) is an insulative, ablative TPS material that is easily 
machined to custom shapes. Machined SIRCA tiles are then applied directly to the spacecraft. The nomenclature of 
SIRCA-15F refers to SIRCA tile manufactured from 0.19 g/cm 3 (12 lb/ft 3 ) bulk density fibrous refractory ceramic 
insulation (FRCI-12), with a final SIRCA bulk density of 0.264 ± 0.024g/cm 3 (16.5 ± 1.51b/ft 3 ). 25 For this study, a 
virgin density of 0.259 g/cm 3 (16.2 lb/ft 3 ) was used for modeling material response. SIRCA has been used 
successfully as an aftbody TPS material for Mars Pathfinder 26 and Mars Exploration Rovers. 27 For these probes, the 
Backshell Interface Plate (BIP) was covered with SIRCA tiles. 

The material stackup for this MER consists only of SIRCA, and the FIAT SIRCA model is version 1.00. The 
SIRCA MER is given in Eq. 4, and it has a standard deviation of 7.4% to FIAT prediction. Specifics of the MER are 
given in Table 8. Shown in Fig. 6 is the MER fit of the data, and a scatter plot of Goodness of Fit with FIAT- 
predicted TPS thickness is given in Fig. 7. 


TH = 0.5281 



0.5416 


( 4 ) 


Table 8. SIRCA MER Details 


Variable 

Values 

Maximum allowable CW heat flux, W/cm 2 

100 

Recession, cm 

none 

Accuracy to FIAT at one SD 

7.4% 

Largest under prediction of FIAT (% of FIAT) 

16.6 

Trajectories used for correlation 

835 

Minimum thickness, cm 

0.614 
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Figure 6. SIRCA MER banded by one standard deviation 
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FIAT Predicted SI RCA Thickness, cm 

Figure 7. SIRCA MER Goodness of Fit 

2. Acusil II 

Acusil II is a silicone foam material that can be molded to cover complex surface geometries. 28 This material 
was used on the Mars Science Laboratory 29 aftbody to protect the parachute cone cover (PCC) and backshell 
interface plate (BIP). It has a nominal density of about 0.25 g/cm 3 (16 lb/ft 3 ). 

The Acusil II MER is given in Eq. 5, and it has a standard deviation of 7.6% to FIAT prediction. Specifics of the 
MER are given in Table 9. Shown in Fig. 8 is the MER fit of the data, and a scatter plot of Goodness of Fit with 
FIAT-predicted TPS thickness is given in Fig. 9. 



TH = 0.623 


(HL ^ 0S697 


( 5 ) 
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Table 9. Acusil II MER Details 


Variable 

Values 

Maximum allowable CW heat flux, W/cm 2 

100 

Recession, cm 

none 

Accuracy to FIAT at one SD 

7.6% 

Largest under prediction of FIAT (% of FIAT) 

15.1 

Trajectories used for correlation 

835 

Minimum thickness, cm 

0.614 



Figure 8. Acusil II MER 



Figure 9. Acusil II MER GoF 

3. SLA-561V 

SLA-561V is a filled cork silicone packed in a phenolic fiberglass honeycomb. 30 It contains a room temperature 
vulcanizing (RTV) silicone resin, granulated cork, silica and phenolic microballoons, and chopped Refrasil fibers. 
The nominal density of the virgin material is about 0.256 g/cm 3 (16 lb/ft 3 ). 
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For this MER, SLA-561V is listed as an aftbody heat shield material because forebody peak heating is at least 
151 W/cm 2 (cold-wall). However, this material has been used on the forebody TPS for all Mars probes with the 
exception of Mars Science Laboratory. 31 

The SLA-561V MER is given in Eq. 6, and it has a standard deviation of 8.5% to FIAT prediction. Specifics of 
the MER are given in Table 10. Shown in Fig. 10 is the MER fit of the data, and a scatter plot of Goodness of Fit 
with FIAT-predicted TPS thickness is given in Fig. 11. 

TH = 0.0064 (|^) 2 + 0.0961 + 0.3322 (6) 


Table 10 SLA-561V MER Details 


Variable 

Values 

Maximum allowable CW heat flux, W/cm 2 

100 

Recession, cm 

none 

Accuracy to FIAT at one SD 

8.5% 

Largest under prediction of FIAT (% of FIAT) 

15.7 

Trajectories used for correlation 

835 

Minimum thickness, cm 

0.454 



Figure 10. SLA-561V MER 
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Figure 11. SLA-561V MER Goodness of Fit 

4. LI-900 

LI-900 tiles are made from 99.9% pure silica glass fibers 32 and are 94% air by volume. The nominal density of 
LI-900 is 0.144 g/cm 3 (9 lb/ft 3 ). These tiles were used on the lower surface of the Space Shuttle Orbiter. 33 For this 
MER, the LI-900 surface is uncoated. Having no surface coating results in a conservative sizing estimate due to an 
uncoated surface having a lower surface emittance 34 and consequently lower radiation from the vehicle surface 
during entry. There was no radiation to the vehicle aftbody, so absorptivity did not have any effect on surface 
heating. 

Although the ratio of heat load and the square of entry velocity is kept as an independent variable, the LI-900 
MER is divided by peak heat flux into two separate equations. This was done to increase MER accuracy. These 
relations are given in Eq. 7. The MER has a standard deviation of 14.0% to FIAT prediction. Specifics of the MER 
are given in Table 7. Shown in Fig. 12 is the MER fit of the data, and a scatter plot of Goodness of Fit with FIAT- 
predicted TPS thickness is given in Fig. 13. 

TH = 0.6961 (-) q cw < 10 

TH = -0.0306 + 0.5896 + 0.6739 q cw > 10 

Table 1 LI-900 MER Details 


Variable 

Values 

Maximum allowable CW heat flux, W/cm 2 

75 

Recession, cm 

none 

Accuracy to FIAT at one SD 

14% 

Largest under prediction of FIAT (% of FIAT) 

18.7 

Trajectories used for correlation 

798 

Minimum thickness, cm 

0.686 


( 7 ) 
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a) qmax <= lOW/cm 

Figure 12. LI-900 MER 
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Figure 13. LI-900 MER Goodness of Fit 
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V. PICA Arcjet Testing Database 

For all the trajectories considered in this work, the highest predicted heating is over 3700 W/cm 2 . PICA has 
never been tested under such high heating conditions (see Fig. 14), and it is unlikely to survive such a harsh 
environment. Figure 14 shows the heat flux for several successful PICA arc jet tests 35,36,37 . The PICA MER is 
limited to a cold-wall peak heating of 1200 W/cm 2 , and M-SAPE has been configured so that a user is notified if a 
trajectory exceeds its MER limit. As a comparative note, it was estimated that a peak heat flux of about 1220 W/cm 2 
was experienced on the PICA heat shield of the Stardust 38 capsule. 
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Figure 14. Successful Stagnation Tests of PICA in which no spallation was observed. Heat flux values are 
cold-wall. 


VI. Conclusion 

Mass estimating relationships have been presented for the vehicle forebody ablators PICA and Carbon Phenolic 
atop ACC, and for the backshell materials SIRCA, Acusil II, SLA-561V, and LI-900. These MERs are accurate to 
FIAT prediction between 7% to 14% at one standard deviation. The most any MER can under predict FIAT is 
18.7% Applications include quick estimates of TPS mass during early stages of vehicle design. These MERs have 
been integrated into M-SAPE and used with FIAT as an initial estimate of required material thickness to speed up 
sizing estimates. When using these MERs, care needs to be taken so that sizing environments, such as peak heating, 
are within the capabilities of the material. Using these MERs has reduced design turnaround times for a possible 
Earth entry configuration from weeks to hours. 
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